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Abstract 

I introduce a gauge invariant decomposition of the nucleon spin into quark 
helicity, quark orbital, and gluon contributions. The total quark (and hence 
the quark orbital) contribution is shown to be measurable through virtual 
Compton scattering in a special kinematic region where single quark scatter- 
ing dominates. This deeply-virtual Compton scattering (DVCS) has much 
potential to unravel the quark and gluon structure of the nucleon. 
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The spin structure of the nucleon reflects interesting non-perturbative physics in Quan- 
tum Chromodynamics (QCD). ^From the recent data on polarized deep- inelastic scattering 

one finds that about 20 ± 15% of the nucleon spin is carried by quark spin or helicity [0. 
Natural questions are then where is the remainder of the nucleon spin? How to measure or 
calculate it? This Letter attempts to provide an answer to them. 

Intuitively, the candidates for the "missing" spin are the quark and gluon orbital angular 
momenta and gluon helicity. In QCD, they can be identified with matrix elements of certain 
quark-gluon operators in the nucleon state || . The problem, however, is that these operators 
take free-field expressions and are not gauge-invariant in an interacting gauge theory. Hence 
it is doubtful that their matrix elements have any experimental significance, although they 
can be calculated in theory, for instance on a lattice, with a fixed gauge. 

In this Letter I show that there exists a gauge invariant decomposition of the QCD 
angular momentum operator into quark and gluon contributions. The quark part can be 
separated further into the usual quark helicity plus the gauge-invariant orbital contribution. 
There exists, however, no gauge-invariant separation of the gluon part into helicity and 
orbital contributions, although high-energy scattering favors such a separation in the light- 
like gauge and infinite momentum frame. The gauge-invariant quark and gluon contributions 
to the nucleon spin is shown to asymptotically approach ratio 16 : 3nf, where nf is the 
number of active fermion flavors. This result is incidentally the same as what Hoodbhoy, 
Tang and this author have derived in a gauge non- invariant formulation 0]. The gauge- 
invariant expression for the angular momentum operator allows one to calculate meaningfully 
fractions of the nucleon spin carried by quarks and gluons. Furthermore, it allows them to 
be measured in deeply-virtual Compton scattering (DVCS) in which the virtual photon 
momentum approaches the Bjorken limit. DVCS gives an access to a new class of nucleon 
observables — the off-forward parton distributions — which are a generalization of ordinary 
parton distributions and elastic form factors. Therefore, effectively DVCS provides a new 
ground to explore the quark and gluon structure of the nucleon. 

The angular momentum operator in QCD is defined according to the generators of 
Lorentz transformation, 

ji = l _ t nk r d 3 xM ojk ^ ^ 



where M ^ is the angular momentum density, expressible in terms of the energy-momentum 
tensor T^ v through 

T^ U has the Belinfante-improved form and is symmetric, gauge- invariant, and conserved ||. 
It can be separated into gauge-invariant quark and gluon contributions, 



where the quark part is, 



and the gluon part is, 



Tf = -[ip^iD^ip + ^iS^V] , (4) 
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= X -y v F 2 - F^ a F u a , (5) 
y 4 

where (fiv) denotes symmetrization with respect to //, v indices. I will ignore the issues of 
gauge fixing and trace anomaly ||, as they do not affect the following discussion. 

For the above equations, one sees that J can be written as a gauge invariant sum, 
</qcd = J g + J g , where 

JL = h lJk I d 3 x(T™x j - T°ix k ) . (6) 



Q,9 2 J q ' 9 q,g 

In pure gauge theory, J g by itself is a conserved angular momentum charge, generating spin 
quantum numbers for glueballs. It is clear that J q and J g are interaction-dependent and 
thus differ from the corresponding expressions in free-field theory. 

To understand the physical content of the gauge invariant J q and J g , one can re-express 
them using QCD equations of motion and superpotentials ||||. After some algebra, one 
finds, 



J q = I d 3 x ^[775 + (x x iD)]if} , 

J g = f d 3 x (xx(Ex B)) , (7) 



where color indices are implicit. One could have guessed the result with going through the 
formal derivation: The quark angular momentum contains the usual quark helicity operator 
plus a gauge-invariant orbital contribution. The gluon angular momentum is constructed 
from the Poynting momentum density E x B. It is a bit surprising though that the simple 
form works for an interacting gauge theory. 

At this point, it is instructive to compare the above gauge invariant form of the angular 
momentum operator with the free- field expression 0. Because the interaction between 
quarks and gluons contains no derivative, one can write Jqcd = J' q + J' g , where J' are 
interaction- independent (except in the definition of E), 

J' q = J d 3 x ip^ 775 + (x x id) ip , 

J' g = J d 3 x [(E xA)- Ei(x x d)Ai] . (8) 

Here each term has a straightforward interpretation: Both the quark and gluon angular 
momenta are sum of spin and orbital contributions. In presence of gauge interactions, J' q and 

J g , as well as the individual terms in them except the quark helicity, are gauge- dependent. 
For this reason, it is difficult to find experiments to measure the other contributions to the 
nucleon spin. 

The gluon angular momentum J g does not admit further gauge-invariant decomposition 
as spin and orbital contributions, contrary to the quark case. This is because the spin space 
of gluons coincides with the ordinary space and time, and therefore under spatial rotation, 
the orbital and the spin representations of the Lorentz group mix. If one disregards the 
issue of gauge invariance, J g can be viewed as a sum of three terms: two terms in J' g in 
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Eq. (8) and an interaction-dependent term —g J d 3 ip^x x Aip. The term S g = J d 3 xE x A 
has the simple interpretation as the spin of gluons in the A = gauge 0,[7|,||]. in t ne finite 
momentum frame (IFM), the gauge condition becomes A + = and the rotational generators 
are constructed from the density M +tJ '. The nucleon matrix element of S g in the light-like 
gauge and IFM is measurable in high-energy scattering. In fact, the first moment of the 
polarized gluon distribution Ag(x) gives pL0], 



Ag{x)dx 2S= (PS\6\PS) 



(9) 



where O is a gauge-invariant operator which reduces to S g in the A + = gauge and IFM . 
To maximally utilize this piece of experimental information, one can define the gluon orbital 
angular momentum L g as the difference J g — S g . The nucleon matrix element of L g in the 
light-like gauge and infinite momentum frame can be deduced from the matrix elements of 
S and J g and might offer some insights on the spin structure of the nucleon. 

According to Eq. (|(J), the Q 2 evolution of the quark and gluon contributions to the 
nucleon spin is the same as that of matrix elements of quark and gluon operators. The 
reason is quite simple: forming spatial moments of T£ v and T^ v does not change the short- 
distance singularity of the operators. On the other hand, in Ref. it was shown that 
the matrix elements of J' q and J' g have the same leading-log evolution as the quark and 
gluon contributions to the nucleon momentum. If both results above are consistent, the 
interaction-dependent term, —g J d 3 i/;'x x Aip, shall not affect the leading-log evolution in 
the light-like gauge. Indeed, an explicit calculation confirms this. If one defines, 



J q , g {Q 2 ) 2S= (PS\J q , g (Q 2 )\PS) 



(10) 



the leading-log Q 2 dependence is simply, 



3n 



' lnQg/A 2 * 
In Q 2 /A 2 / 

' lnQg/A 2 X 
216 + 3n/ ' VlnQ 2 /A 2 



2 16 + 3n f 
1 16 



+ 



+ 



2(16+3n / )/(33-2n / ) 



2(16+3n / )/(33-2n f ) 



UQl) 



UQl) - 



1 3n 



f 



2 16 + 3n f _ 

1 16 

216 + 3n f 



where n/ is the number of quark flavors (J q + J g — 1/2). As Q 2 ^ oo, the partition of 
the nucleon spin between quarks and gluons approaches the ratio 16 : 3rif, the same as the 
asymptotic partition of the nucleon momentum derived by Gross and Wilczek 0. 

The gauge-invariant form of the QCD angular momentum operator allows one to mean- 
ingfully calculate and measure the fraction of the nucleon spin carried by quarks and gluons. 
Recently, Balitsky and I have estimated J q ^ g at the scale of 1 GeV using the QCD sum rule 
method |10| |. To see how they can be measured in an experiment, I define the form factors 
of the quark and gluon energy-momentum tensors, 



[P'\T^ g \P) = f/(P')[A,s( A2 )7 (Atpl/) +B q , g (A 2 )P^ia^ a A a /2M 

+0, ig (A 2 )(AW - g^A 2 )/M + C qtg {A 2 )g^M]u{P) 



(12) 



where P^ = (P M + P M ')/2, A M = P^' — P M , and U(P) is the nucleon spinor. [An analogous 
equation for the full tensor was considered by Jaffe and Manohar H].] The barred form factor 
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arises from non-conservation of the tensor currents. Taking the forward limit for /x = and 
integrating over 3-space, one finds that A qtg (0) give the momentum fractions of the nucleon 
carried by quarks and gluons (A q (Q) + A g (0) = 1). On the other hand, substituting the 
above into the nucleon matrix element of Eq. (|^), one finds, 

J q , 9 = I [4r*(°) + ■ (13) 

Thus, to find the quark and gluon contributions to the nucleon spin, one has to measure the 
P-form factor, which is analogous to the Pauli form factor for the vector current. 

Since there is no fundamental probe that couples to the quark and gluon energy- 
momentum tensors (the graviton does, but only to the sum), it appears hopeless to measure 
the form factors. On the other hand, T^ v does appear in the operator product expansion 
(OPE) for the product of vector currents TJ a (£)Jp(0), and thus B qtg (0) is accessible through 
deep- inelastic sum rules, like what A q)g (0 is measured. The complication, however, is that 
the B qj g(Q) term does not contribute to the forward matrix element and so the usual inclu- 
sive deep-inelastic process is useless. A way to get around this is to measure the off-forward 
matrix element of TJ a (£)Jp(Q) and extrapolating the form factors to the forward limit. The 
natural process to do this is Compton scattering |TT|. To ensure there is an OPE, one has to 
have one of the photons far off-shell and single-quark scattering dominating the process. To 
emphasize this special kinematic region, I call the process deeply- virtual Compton scatter- 
ing (DVCS). In the remainder of this Letter, I summarize the main features of this process 
without going through any detailed proof. 

Consider virtual Compton scattering with a virtual photon of momentum absorbed by 
a nucleon of momentum P M , and an outgoing real photon of momentum q ,fl = — A M , and a 
recoil nucleon of momentum P' M = P M + A M . The deeply- virtual kinematics refer to q^ in the 
Bjorken limit, namely, Q 2 = —q 2 — > oo, P ■ q — > oo, and Q 2 /P ■ q finite. To ensure a reliable 
extrapolation to A M = 0, the components of A M shall be as small as possible. However, as 
will be clear below, the best one can do is to constrain them on the order of the nucleon mass. 
Given the kinematics above, it is easy to show that the only dominant scattering subprocess 
involves a single quark absorbing the deeply-virtual photon and subsequently radiating a 
real photon and falling back to the nucleon. Other subprocesses are down by at least a factor 
of 1/Q 2 and are negligible. I henceforth concentrate on the dominant subprocess shown in 
Fig. 1. 

To calculate the scattering amplitude, it is convenient to define a special system of 
coordinates. I choose g M and P M = (P + P') M /2 to be collinear and in the z direction. 
Introduce two light-like vectors, p^ = A(l, 0,0,1) and = (1, 0, 0, — 1)/(2A), with p 2 = 
n 2 = 0, p ■ n = 1, and A arbitrary. I expand other vectors according to n fl and transverse 
vectors, 

P»=p» + (M 2 /2)ri l 
q» = + (Q 2 /2£K 

A^ = -£(p M - M 2 /2n tM ) + A^ (14) 

where M 2 = M 2 - A 2 /4 and f = Q 2 /(2P ■ q). The Compton amplitude = 
i J d i ye~ iq ' y (P'\TJ v {y)J^{Q)\P), where v and /i are the polarization indices of the initial 
and final photons, is, 



5 



T^(P,g,A)= _I(p<V+pV-(r) / rfa; f— 



+ X + £/2 + ie, 



X 



A 2 , A • n)U{P')iU{P) + P(x, A 2 , A • n)U(P' 



2M 



U{P) 



X 



A • 71, 



#(ar, A 2 , A • n)^(P') ^ 75 ^(P) + E(x, A 2 , A • n)— — U(P') l5 U(P) 



2M 



• (15) 



Thus only 4 of the 12 helicity amplitudes survive the Bjorken limit ||12|| . All photon helicity- 
flipping and longitudinal photon amplitudes vanish. H, H, E and E are new, off-forward, 
twist-two parton distributions defined through the following light-cone correlation functions, 

^e iXx (P'\i>(-Xn/2)YH^/2)\P) = H(x, A 2 , A • n)U(P , )-fU{P) 
2tc 

+E(x, A 2 , A ■ n)U(P')^j^U(P) + ... 
^e iXx {P'\${-\n/2)ri^{Xn/2)\P) = H{x, A 2 , A • n)U{P')^ lh U{P) 

Z7T 

+E(x, A 2 , A ■ nmP'Y-^UiP) + ... (16) 

where I have negelected the gauge link and the dots denote higher-twist distributions. ^From 
the definition, H and H are nucleon helicity-conserving amplitudes and E and E are helicity- 
flipping. 

The off-forward parton distributions have the characters of both ordinary parton distri- 
butions and nucleon form factors. In fact in the limit of A^ — > 0, we have 



H(x,0,0) = f 1 (x), H{x,0,0)= 9l (x) 



(17) 



where f± (x) and g\ (x) are quark and quark helicity distributions. On the other hand, forming 
the first moment of the new distributions, one gets the following sum rules, 



J dxH(x. 


A 2 . 


A 


n) = F l {A 2 ) 


J dxE(x, 


A 2 . 


A 


n) = P 2 (A 2 ) 


J dxH{x. 


A 2 . 


A 


n) = G A {A 2 ) 


JdxE{x. 


A 2 . 


A 


n) = G P (A 2 ) 



(18) 

where Pi and F 2 are the Dirac and Pauli form factors and Ga and Gp are the axial- vector 
and pseudo-scalar form factor. The most interesting sum rule relevant to the nucleon spin 
is, 



J dxx[H(x, A 2 , A • n) + E(x, A 2 , A • n)} = A q (A 2 ) + B q (A 2 ) 



(19) 
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where luckily the A • n dependence, or C q (A 2 ) contamination, drops out. Extrapolating the 
sum rule to A 2 = 0, the total quark (and hence quark orbital) contribution to the nucleon 
spin is obtained. By forming still higher moments, one gets form factors of various high-spin 
operators. 

It is important to comment on practical aspects of the experiment. First of all, from the 
cross section, one finds that E and H can be measured either in unpolarized scattering, or in 
electron single-spin asymmetry through interference with the Bethe-Heitler amplitude jI2 



or in polarized electron scattering on a transversely polarized target. A detailed examination 
of various possibilities, together with some numerical estimates will be published elsewhere. 
Second, the DVCS cross section is down by an order of a em compared with the deep-inelastic 
cross section, but has the same scaling behavior. So the cross section is measurable, but 
statistics would be a challenging requirement. The ideal accelerator for the experiment is 
ELFE ]nj. Finally, the extrapolation of A 2 from order M 2 to requires dispersive study of 



the form factors of the tensor currents. The B q (A ) form factor is dominated by resonances 
in the exotic 1 h channel. 

The off-forward parton distributions can be defined for quark helicity-flip (chiral-odd) 
correlations, for higher twists, and for gluons. DVCS provides one process to access to these 
distributions. There are other processes one can consider to measure them. For instance, the 
diffractive p or J/ip production studied recently by Brodsky et al. 0] can be used to measure 
the off-forward gluon distributions. Thus there is now a new territory to explore the quark 
and gluon structure of the nucleon besides the traditional inclusive (parton distributions) 
and exclusive (form factors) processes. 
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P P' = P - A 

FIG. 1 . Dominant scattering process in deeply- virtual Compton scattering 



